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Recently, a measure for the non-Markovian behavior of quantum processes in open systems has 
been developed which is based on the quantification of the flow of information between the open 
system and its environment [Phys. Rev. Lett. 103, 210401 (2009)]. The information flow is connected 
to the rate of change of the trace distance between quantum states which can be interpreted in terms 
of the distinguishability of these states. Here, we elaborate the mathematical details of this theory, 
present applications to specific physical models, and discuss further theoretical and experimental 
implications, as well as relations to alternative approaches proposed recently. 

PACS numbers: 03.65.Yz, 03.65.Ta, 42.50.Lc 



I. INTRODUCTION 

A Markov process in the evolution of an open quantum 
system typically gives rise to a quantum dynamical semi- 
group for which the most general representation can be 
written in the Lindblad form p] |2] . There exist however 
complex systems for which this relatively simple descrip- 
tion of the open system dynamics in terms of a Markovian 
master equation fails to give a comprehensive picture of 
the dynamics [3] . Thus in many realistic physical systems 
the Markovian approximation of the dynamics gives an 
overly simplified picture of the open system evolution and 
a more rigorous treatment of the dynamics is required. 

To give insights into the nature of non-Markovian ef- 
fects many analytical methods and numerical simulation 
techniques have been developed in recent years (see, for 
example, Refs. [HH7|). Non-Markovianity manifests it- 
self in the different approaches in a variety of ways and 
there exists no general recipe for comparing the degree of 
non-Markovianity in different physical models. In order 
to give a general quantity determining the degree of non- 
Markovian behavior in the open system dynamics, one 
has to rigorously define what makes a dynamical map 
non-Markovian. 

Here, we discuss a recently proposed measure for the 
degree of non-Markovian behavior which is based on the 
trace distance between quantum states [TH]. The trace 
distance describes the probability of successfully distin- 
guishing two quantum states and the change in the trace 
distance of two open system states can be interpreted as 
a flow of information between the system and the envi- 
ronment. When the trace distance decreases information 
flows from the system into the environment, while an in- 
crease of the trace distance signifies a backflow of infor- 
mation from the environment to the system. Markovian 



processes tend to continuously decrease the distinguisha- 
bility between any two states of the open system, i.e., in- 
formation flows continuously from the system to the envi- 
ronment. The condition which defines a non-Markovian 
dynamical map is that the map allows an information 
flow from the environment to the system and, therefore, 
allows the system to gain information about its former 
state. This condition for a non-Markovian map leads to 
a rigorous and general definition of a measure for the 
degree of non-Markovianity in open quantum systems. 

In Sec. |H] we construct the measure for non- 
Markovianity and discuss its properties for some general 
classes of quantum processes in open systems. It is shown 
that the non-divisibility of the dynamical map is neces- 
sary for the process to be non-Markovian. Hence, the 
measure vanishes for quantum dynamical semigroups and 
for time-dependent Markov processes. We also demon- 
strate that the appearance of negative rates in the quan- 
tum master equation is a necessary condition for non- 



Markovianity. In Sec. Ill we illustrate the determination 
of the measure for a two-level system and for a A-type 
atom in a cavity. Section |IV| contains a detailed discus- 
sion of several alternative ways for defining a measure 
for non-Markovianity. Moreover, we present possible ex- 
perimental strategies for the detection of non-Markovian 
effects. The conclusions are drawn in Sec. IV] 



II. THE MEASURE FOR 
NON-MARKOVIANITY 

A. Construction of the measure 

To construct the measure for non-Markovianity we 
need a measure for the distance between any pair of quan- 
tum states represented by density matrices p\ and pi. 
Such a measure is given by the trace distance, which is 
defined as 
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where the modulus of an operator A is defined by \A\ — 
\A4t A The trace distance D yields a natural metric on 
the state space and satisfies < D < 1. It has many nice 
properties that make it a useful measure for the distance 
between quantum states p^. First, the trace distance is 
preserved under unitary transformations U, 



D(U Pl U\U P2 U^) = D{p uP2 ). 



(2) 



Second, all completely positive and trace preserving 
(CPT) maps <E> (trace preserving quantum operations) 
are contractions for this metric, 



D(<f>px,<f>p 2 ) < D{p uP2 ). 



(3) 



Third, the trace distance has a physical interpretation 
as a measure of state distinguishability. Suppose Alice 
prepares a quantum system in the state p\ with proba- 
bility 1/2, and in the state pi with probability 1/2. She 
gives the system to Bob, who performs a measurement 
to distinguish the two states. The maximal probability 
that Bob can identify the state given to him is [2"U] 



Pmax = 2 I 1 + D iPUP2)] ■ 



(4) 



Hence, the trace distance represents the maximal bias in 
favor of the correct state identification which Bob can 
achieve through an optimal strategy. For example, if 
pi and pi have orthogonal supports the trace distance 
becomes D(p\,p2) = 1 and thus p max = 1, which means 
that Bob is able to distinguish the states with certainty. 

The change in the distinguishability of states of an 
open system can be interpreted as a flow of information 
between the system and the environment. We consider 
here quantum processes given by a dynamical CPT map 
<I>(i, 0) which transforms the initial states p(0) at time 
zero to the states p(t) at time t > 0, 



p(0)H>p(t) = $(t,0)p(0). 



(5) 



When such a quantum process reduces the distinguisha- 
bility of states, information is flowing from the system 
to the environment. Likewise, the increase of the dis- 
tinguishability signifies that information flows from the 
environment to the system. The invariance under unitary 
transformations ^ indicates that information is pre- 
served under the dynamics of closed systems. The con- 
traction property of Eq. ^ guarantees that the maximal 
amount of information the system can recover from the 
environment is the amount of information earlier flowed 
out the system. 

The basic idea underlying our construction for the 
measure of non-Markovianity in a quantum process is 
that for Markovian processes information flows contin- 
uously from the system to the environment. In order 
to give rise to non-Markovian effects there must be, for 
some interval of time, an information flow from the en- 
vironment back to the system. The information flowing 
from the environment back to the system allows the ear- 
lier states of the system to have an effect on the later 



dynamics of the system, i.e., it allows the emergence of 
memory effects. 

We define the rate of change of the trace distance of a 
pair of states by means of 



<r(t,Pi,2{°)) = j t D(pi{t),p 2 (t)), 



(6) 



where pi,2(t) — 3>(t, 0)pi.2(0). For a non-Markovian pro- 
cess described by a dynamical map <&(£,()), information 
must flow from the environment to the system for some 
interval of time and thus we must have a > for this 
time interval. A measure of non-Markovianity should 
measure the total increase of distinguishability over the 
whole time evolution, i.e., the total amount of informa- 
tion flowing from the environment back to the system. 
This suggests defining the measure 7V($) for the non- 
Markovianity of the quantum process 0) through 



J\f(§) = max / 

Pi,2(0) J CT>0 



dta(t,p lt2 (0)). 



(7) 



The time integration is extended over all time intervals 
(aj, hi) in which a is positive and the maximum is taken 
over all pairs of initial states. Due to Eq. ^ the measure 
can be written as 



7V($)= max Y J [ D {piQ>i),p2{b i ))-D(p 1 {a i ),p 2 {a i ))\. 

Pi, 2(0) ^ 

% 

To calculate this quantity one first determines for any 
pair of initial states the total growth of the trace dis- 
tance over each time interval (a,,^) and sums up the 
contribution of all intervals. Af($) is then obtained by 
determining the maximum over all pairs of initial states. 
While it may be difficult to derive an analytical expres- 
sion for the measure defined in Eq. ([8]), the numerical 
evaluation of the measure is relatively easy provided the 
dynamical map is known explicitly. We will discuss in 
Sec. |III| the determination of Af($) for some specific ex- 
amples. 



B. Classification of quantum processes 

Having defined our measure for non-Markovianity we 
discuss in this section the properties of this measure for 
some general classes of quantum processes. Specific phys- 



ical systems will be investigated in Sec. Ill 



1. Divisible maps 

The dynamical map Q(t, 0) is defined to be divisible if 
for all t, t > the CPT map $(< + r, 0) can be written as 
composition of the two CPT maps $(t + r, t) and <&(i,0), 



$(< + t,0) = *(* + T, <)*(*, 0). 



(9) 



We note that this definition differs slightly from the usual 
definition of divisibility according to which a CPT map 
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A (quantum channel) is said to be divisible if there exist 
CPT maps Ai and A2 such that A = A1A2, where it is as- 
sumed that neither Ai nor A 2 is a unitary transformation 
[2"T] . In Eq. ([9]) the left-hand side as well as the second 
factor on the right-hand side are fixed by the given dy- 
namical map. Hence, Eq. ([9| requires the existence of a 
certain linear transformation $(i + t, t) which maps the 
states at time t to the states at time t + t and represents 
a CPT map (that may be a unitary transformation) for 
all t and all r. There are many quantum processes which 
are not divisible. For instance, if <&(i,0) is not invertible, 
a linear map $(t + r, t) which fulfills Eq. Q may not 
exist. Moreover, even if a linear map Q(t + r,t) satisfy- 
ing Eq. ([£]) does exist, this map needs not be completely 
positive, and not even positive. 

We claim that all divisible dynamical maps are Marko- 
vian. To prove this statement suppose that $(i, 0) is di- 
visible. For any pair of initial states /3i,2(0) we then have 



Pi, 2 (i + r) - $(t + r,iW 2 (t). 



(10) 



Since $(< + r, t) is a CPT map we can apply the contrac- 
tion property ^ to obtain: 



D( Pl (t 



,P2(t + r))<D( Pl (t),p 2 (t)). 



(11) 



This shows that for all divisible dynamical maps the trace 
distance decreases monotonically, i.e., a(t, pi,2(0)) < 
and, therefore, A/"($) = 0. Thus, we conclude that all di- 
visible processes are Markovian and that non-Markovian 
processes must necessarily be described by a nondivisible 
dynamical map. 



2. Quantum dynamical semigroups 

The prototype of a Markovian dynamics is provided by 
a Markovian master equation for the density matrix, 



dt 



P (t) = Cp(t), 



with a generator in Lindblad form [TJ [2] 



Cp 



-i\H, 
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(12) 



, (13) 



involving a time-independent Hamiltonian H as well as 
time-independent Lindblad operators Ai and positive de- 
cay rates 7^ > 0. Such a master equation leads to a 
dynamical semigroup of CPT maps, $(i, 0) = exp(Ct). 
With $(£ +T, i) — exp(£r) the divisibility condition Q 
is trivially satisfied. Hence, we have AA(4>) = for all 
dynamical semigroups, i.e., for all processes described by 
a master equation in the Lindblad form. 



3. Time- dependent Markov processes 



The divisibility property holds for a much larger class 
of quantum processes than those described by a master 



equation of the form (12). Suppose we have a time-local 



master equation of the form 



Jt p{t)=lC{t)p{t) 



(14) 



with a time-dependent generator IC(t). It can be shown 
that in order to preserve the Hermiticity and trace of the 
density matrix this generator must be of the form [TJ |S] 



K{t)p 



+ 



-i[H(t),p] (15) 
£ 7 i(i) \Mt) P A\(t) - \ {Al(t)Mt),p} 



By contrast to the assumptions in Eq. ( 13 ) the Hamilto- 



nian H(t), the Lindblad operators Ai(t) and the decay 
rates "fi(t) may now depend on time. If the decay rates 
are positive functions, ^(t) > 0, the generator (15) is in 
Lindblad form ( 13 1 for each fixed t > 0. Such a process 



with 7i(t) > may be called time-dependent Markovian 
although the corresponding dynamical map 



$(i,0) = Texp 



dt'fC(t') 



(16) 



does not yield a dynamical semigroup (T denotes the 
chronological time-ordering operator). However, one can 
easily see that the divisibility condition ^ still holds 
because the map 



$(t + r,t) = Texp 



t+T 



db'K,(t') 



(17) 



is CPT for 7i(t) > 0. Thus we can conclude that for 
all time-dependent Markovian processes we again have 
Af(&) = 0. 

We have just seen that a quantum process given by 
the time-local master equation ( 14 ) with positive rates 
leads to a divisible dynamical map. Under certain con- 
ditions the converse of this statement is also true. More 
precisely, if the dynamical map $(i, 0) is divisible with a 
unique map $(t+r, t) depending smoothly on r, then the 
corresponding density matrix p(t) obeys a master equa- 
tion of the form ( 14 ) with positive rates in the generator 
(fl5|). In fact, using p(t + r) = $(£ + t, t)p(t) we find 



d , . 
dt^ 



d 

Jr 



(18) 



and, hence, we obtain the master equation ( 14 ) , where 
the generator is given by 



Kit) 



d 

Jr 



(19) 



Since <!>(£ + r, t) is CPT and satisfies $(t,t) = I, this 
generator must be in Lindblad form for each fixed t, i.e., 
it must have the form (15) with ji(t) > 0. 
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4- Non-Markovian processes 

The measure for quantum non-Markovianity does not 
depend on any specific mathematical representation of 
the dynamics. There are many different such representa- 
tions, e.g., through generalized master equations involv- 
ing a certain memory kernel. However, quantum master 
equations with the time-local structure given by Eqs. ( f 4 ) 



and ( f 5 ) are also very useful for the description of non- 



Markovian processes. It follows from the preceding re- 
sults that in order for such a master equation to yield 
a nonzero measure, 7V($) > 0, at least one of the rates 
7i(i) must take on negative values for some interval of 
time. We emphasize that temporarily negative rates in 
the master equation do in general not lead to a violation 
of the complete positivity of the dynamical map. Many 
examples for time-local master equations with negative 
rates are known in the literature. Further examples will 
be discussed in the next section. 



III. EXAMPLES 
A. Two-level system 

We study the dynamics of a two-level atom with ex- 
cited state |+) and ground state |— ) which is coupled 
to a reservoir of held modes initially in the vacuum 
state. In Ref. [TS] we have described the detuned Jaynes- 
Cummings model, while here we treat the resonant case. 
We will show that the pair of states maximizing the mea- 
sure for non-Markovinity is different in the two cases. 
This demonstrates that the change in both the popula- 
tions and the coherences plays a crucial role in the flow 
of information between the system and the environment. 

The two-level atom model can easily be solved exactly 
[3] and leads to a dynamical map Q(t,0) which can be 
represented in terms of the elements p±±(t) = (±|p(i)|±) 
of the density matrix p(t) as follows, 

p ++ (t) = |G(i)|V+(fj), 

P— (*) = P— (0) + (l-|G(t)| 2 ) P++ (0), 

p+-(t) = G(t)p + _(0), 

p-+(t) = GC(t)p_+(0). (20) 

Here, the function G(t) is defined as the solution of the 
integrodifferential equation 



dt 



G(t) 



dtJit-t^Gih) 



(21) 



corresponding to the initial condition G(0) = 1, where 
f(t — tx) denotes the two-point reservoir correlation func- 
tion (Fourier transform of the spectral density). The map 
(20) is completely positive if and only if \G(t)\ < 1. One 



Eq. ([9]), where the map $(i + t, t) is given by 



P++(t + r) 



G{t + r) 



G(t) 

p—(t + r) = p—(t)+ 1 1- 



P++(t), 



G(t + t) 



G(t) 



/>++(*)> 



p+-{t + r) 
p-+(t + r) 



G(t + r) 
G{t) 

G*(t + T) 

G*{t) 



*>+-(*)> 
P-+(t). 



(22) 



It follows from these equations that a necessary and suf- 
ficient condition for the complete positivity of $(t + t, t) 
is given by \G(t + r)| < |G(t)|. Thus we see that the dy- 
namical map of the model is divisible if and only if \G{t) \ 
is a monotonically decreasing function of time. Note that 
this statement holds true also for the case that G(t) van- 
ishes at some finite time. 

With the help of the above results one can easily derive 
an analytical formula for the time derivative of the trace 
distance, 



<K*,pi,2(o)) 



2|G(i)| 2 a 2 + |6| 2 d 



\b\ 2 dt 



\G(t)\, (23) 



where a = p± + (0) — pJ + (0) denotes the difference of the 
populations and b = p\~ (0) — p\~ (0) the difference of 
the coherences of the initial states. This relation shows 
that the trace distance increases at some point if and 
only if \G(t)\ increases at this point. We conclude that 
the measure for non-Markovianity is positive, A/"($) > 0, 
if and only if the dynamical map is nondivisible. 

A positive measure for non-Markovianity is not only 
linked to a breakdown of the divisibility of the dynamical 
map, but also to the emergence of a negative rate in the 
corresponding master equation (14). In fact, as long as 



G(t) ^ one can write an exact master equation of this 
form with the generator 



(24) 



+7(t) 

where we use the definitions 
j(t) = -25R 

Writing the rate ^(t) as 
7(*) = - 



a-pa + - - {<7+a-,p} 




S(t) = -23? 



\G(t)\ 




\G(t)\ dt 



(25) 



(26) 



can easily check that $(f, 0) can be decomposed as in 



we see that an increase of \G(t) \ and, hence, a breakdown 
of the divisibility leads to a negative rate in the generator 
of the master equation. Thus we find that for the present 
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FIG. 1: The non-Markovianity A/"($) for the damped Jaynes- 
Cummings model as a function of the coupling strength 70. 
Gray dots: 1000 randomly drawn pairs of pure and mixed 
initial states. Black circles: The initial pair given by Eq. ( 29 1 
which leads to the maximum in Eq. Q. 



model a nonzero measure for non-Markovianity is equiv- 
alent to the non-divisibility of the dynamical map and to 
the occurrence of a negative rate in the master equation. 

As an example we consider the case of a Lorentzian 
reservoir spectral density which is on resonance with the 
atomic transition frequency and leads to an exponential 
two-point correlation function 



f( T ) = ^7oAe 



-A|t| 



(27) 



where 70 describes the coupling strength and A the spec- 
tral width (damped Jaynes-Cummings model). Solving 
Eq. (21) with this correlation function we find 



G(t) 



-At/2 



cosh 



dt 



A , [dt 

— smh — 
d V 2 



(28) 



where d = \J\ 2 — 270A. We see that for small couplings, 
70 < A/2, the function \G(t)\ decreases monotonically. 
The dynamical map is thus divisible in the weak coupling 
regime, the rate 7(f) is positive, and the measure for non- 
Markovianity vanishes. However, in the strong coupling 
regime, 70 > A/2, the function \G(t)\ starts to oscillate, 
showing a non-monotonic behavior. Consequently, the 
dynamical map is then no longer divisible and A/"($) > 0. 
We note that in the strong coupling regime the rate 7(4) 
diverges at the zeros of G(t). However, the master equa- 
tion can still be used to describe the evolution between 
successive zeros and, therefore, the connection between a 
positive measure and negative rates in the master equa- 
tion remains valid. 

There is thus a threshold 70 = A/2 for the system- 
reservoir coupling below which 7V(<i>) = 0. We find that 
the measure increases monotonically with increasing cou- 
pling for 70 > A/2. This is illustrated in Fig. [j] The 
maximization over the pair of initial states pi^OJin ex- 
pression Q has been performed here by a Monte Carlo 
sampling of pairs of initial states. Our simulations pro- 
vide strong evidence that the maximum is attained for 



the initial states 

Pi(0) = h) (-1 , P 2 (0) = I(|+) + |-»«+|+H). (29) 

In Ref. jTHj we calculated the measure for the detuned 
Jaynes-Cummings model in the weak coupling limit. In 
this example the maximum of the measure was obtained 
for the initial states pi(0) = |— ) (— | and £2(0) = |+) (+|, 
i.e. for the invariant ground state and the excited state. 
The difference in the maximization for the resonant and 
the off-resonant case arises from the fact that the rate 
at which the populations and the coherences initially de- 
cay is much larger for the resonant case. Consequently, 
the growth of the trace distance occurs after the excited 
state population and the coherences have reached the 
value zero. After this point, the increase of the coher- 
ences yields the dominant contribution to the increase 
of the trace distance. Therefore, the maximal growth of 
the trace distance for the resonant case is reached for the 
invariant state and the state with maximal initial coher- 



B. A- model 

The A-model describes a three-level atom with excited 
state j a) and two ground states \b) and |c) interacting 
off-resonantly with a cavity field. This example allows us 
to demonstrate how the measure for non-Markovianity 
operates in a multi-channel case and how there can ex- 
ist simultaneously positive and negative decay rates for 
different channels. The spectral density we use is 



J(w) 



7o 



A 



2-7T (uj cav - uj) 2 + A 2 ' 



(30) 



where w C av is the resonance frequency of the cavity. Fur- 
ther details and the master equation describing the dy- 
namics of the A-type atom are presented in the Ap- 
pendix. The generator of the master equation is of the 
form of Eq. (151 with two Lindblad operators |6) (a| and 
|c) (a|, and two time-dependent decay rates 71 (t) and 
72(i). 

The detunings of the transition frequencies of the A- 
atom from the cavity resonance frequency are denoted 
by Aj = cjj — w cav . When the detuning parameters Ai 
and A 2 are both sufficiently large, the decay rates 71 (t) 
and 72 (t) get temporarily negative values and this gives 
rise to an information flow from the environment to the 
system. On one hand, the two decay rates 71 (t) and 72(0 
have simultaneous negative regions when Ai = A 2 . On 
the other hand, when Ai ^ A 2 , the decay rates can have 
opposite signs. In this case, the co-operative action of the 
other channel reduces the amount of information flowing 
from the environment to the system. The maximum of 
the measure over the initial states is reached when the 
states are chosen to be |o) (a\ and \b) (b\, or \a) (a\ and 
|c) (c|, depending on which of the channels has more infor- 
mation flow from the environment to the system. When 
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Ai and A2 are such that the channel corresponding to 
the decay rate ji(t) (i = 1 or 2) causes more informa- 
tion flow from the environment to the system we get the 
expression 



a(t) = - Ji(t)p aa (t). 



(31) 



The function p aa (t) is specified in the Appendix. Eq. (31 ) 
shows that the A-system is non-Markovian if one of the 
decay rates 71 (t) or ^(t) takes on negative values. The 
maximization over the the pair of initial states is demon- 
strated in Fig. [2j where the measure was again calculated 
numerically from a large sample of initial states. 
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FIG. 2: The non-Markovianity A/"($) for the A-model as a 
function of the detuning Ai for A2/A = 5 and 70/A = 0.01. 
Gray dots: 1000 randomly drawn pairs of initial states. Cir- 
cles: The initial pair pi(0) = a) (a| and P2(0) = |6) (6|. 
Pluses: The initial pair pi(0) = |a) (a| and P2(0) = |c) (c[. 
At Ai = A2 the pair which yields the maximum in Eq. |7|) 
changes from the latter to the former pair of initial states. 



IV. DISCUSSION 

A. Alternative distance measures 

We have based our definition of the measure of non- 
Markovianity on the trace distance ([!]). An alternative 
measure is obtained if one replaces the trace distance by 
the relative entropy 



S(p 1 \\p 2 ) = Tr[pi(logpi - log/02)] 



(32) 



Using this quantity as a measure for the distance be- 
tween quantum states one is led to a similar interpreta- 
tion as before because the relative entropy also decreases 
under CPT maps |22j . There are however some technical 
problems and limitations in the usefulness of the relative 
entropy which arise from the fact that for many pairs 
Pi and pi the relative entropy becomes infinite [23] and 
thus leads to singularities in the definition of the mea- 
sure. This situation can occur even in the simple case of 
a two-state system, demonstrating the problems of the 
relative entropy concept in defining a general measure for 
non-Markovianity. No such problems occur for the trace 



distance which is well-defined and finite for all physical 
states represented by positive trace class operators. 

Another common measure for the distance between 
two states is the Hilbcrt-Schmidt distance 



Dns(pi,P2) = a/Tt [{pi - p 2 ) 2 ]. 



(33) 



For two-dimensional Hilbert spaces the Hilbert-Schmidt 
distance and the trace distance coincide and correspond 
to the Euclidean distance between the Bloch vectors rep- 
resenting the states (up to numerical factors). However, 
the Hilbert-Schmidt distance is not suitable for a defini- 
tion of non-Markovianity since CPT maps are in general 
not contractions for this metric [24]. Thus, the Hilbcrt- 
Schmidt distance does not provide a natural way to define 
the information flow between system and environment. 



B. Experimental issues 

The exact determination of the measure generally re- 
quires solving the complete reduced dynamics which can 
be a difficult task for more complex systems. However, 
any observed growth of the trace distance is a clear sig- 
nature for non-Markovian behavior and leads to a lower 
bound for N {<$>). The measure for non-Markovianity in- 
troduced here could therefore be useful also for the ex- 
perimental detection of non-Markovianity. 

In an experiment one has to perform a state tomogra- 
phy on different ensembles at different times in order to 
decide whether or not the trace distance has increased. 
Such an experiment also allows the validation of theoret- 
ical models or approximation schemes. Consider a theo- 
retical model predicting er(i, /9^ 2 (0)) > for some inter- 
val t G (ti,t 2 ) and for some pair of initial states pi }2 (0). 
In the experiment one should then detect the increase 
of the trace distance between the states pi(t) and p 2 (t) 
in this time interval. This type of experiment could be 
based, e. g., on the recent proposal to use a trapped ion to 
study quantum Brownian motion in the non-Markovian 
regime |25j . The explicit experimental implementation 
of this system can be done, e. g., by using reservoir en- 
gineering techniques |26j or by using the trapped ion as 
a quantum simulator for non-Markovian dynamics [27] . 
One of the possibilities here to detect non-Markovianity 
is to prepare the ion in various Fock states, and to study 
the trace distance dynamics as described above. 

A great advantage of the present approach is that 
it also allows to plan experiments for testing non- 
Markovianity without knowing the properties of the en- 
vironment or the system-environment interaction. The 
interactions and environmental properties can be quite 
difficult to model in an experimental setup. By perform- 
ing a state tomography for two states of the open system 
under study at many different times, one can determine 
whether there has been any increase in the trace distance 
and, hence, non-Markovian behavior in the dynamics. 
From this information one can conclude whether or not 
non-Markovian effects are crucial in the dynamics and 
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in this way also gain some knowledge on the nature of 
the environment and the interactions. An example un- 
der active investigation, where nevertheless a complete 
characterization of the environment is still missing and 
where non-Markovianity could play a role, is given by the 
energy transfer in photosynthctic systems |28j . 



C. Other approaches to non-Markovianity 

Recently, other interesting approaches to the charac- 
terization and quantification of non-Markovianity have 
been proposed. The measure suggested in Ref. [5] quanti- 
fies non-Markovianity in terms of the minimal amount of 
noise required to make a given quantum channel Marko- 
vian. The most important difference to our approach is 
that this measure is based on the properties of the dy- 
namical map at a given time, i.e., on the properties of the 
quantum channel represented by a snapshot of the time 
evolution. Hence, this approach assesses to what extend 
the dynamical map at each fixed time to deviates from 
an element of a Markovian process. The fundamental dif- 
ference between the notion of non-Markovianity used in 
Ref. [5] and ours can be seen from the following simple ex- 
ample. We consider the dynamical map 0) of a two 
state system undergoing a pure de- and re-phasing dy- 
namics which is given by (using the notation of Sec. Ill A ) 



p++(t) = p++(0), p__(t) = p__(0), 

p+-(t) = g(t) P+ -(0), p-+(t) = g(t)p-+(0), (34) 

where the function g[t) = 1 [1 + cos 2 cot] describes a peri- 
odic oscillation of the coherences. The trace distance for 
this model is given by 



D( Pl (t),p 2 (t))= ^ 



>(t)|6|* 



(35) 



where a = pj + (0)-p£ + (0) and b = pj~ (0)-pJ"(0). For 
b 7^ the trace distance thus oscillates periodically and, 
hence, A/"($) = +oo according to the definition ^ of 
our measure. On the other hand, the non-Markovianity 
in the sense of Ref. [5] is zero because for any fixed to the 



dynamical map (34) can be written as an element of a 
Markovian semigroup: $(to,0) = exp(£) with the Lind- 
blad generator Cp = ^(03/303 — p), where V = —h\g(to). 

A further interesting measure proposed recently [5] is 
closely connected to the measure discussed here. In fact, 
the measure of Ref. [5] quantifies deviations from the di- 
visibility of the dynamical map. As we have seen, the 
non-divisibility of the dynamical map is a necessary con- 
dition for to be nonzero. However, we conjecture 
that our notion of non-Markovianity and the one used 
in [6] are not strictly equivalent, i.e., that there are non- 
divisible maps with Af(&) = 0. Further considerations 
concerning this point will be published elsewhere. 



V. CONCLUSIONS 

We have constructed a measure Af($>) for the non- 
Markovianity of quantum processes in open systems in 
terms of the information flowing from the environment 
to the system during the time evolution. The flow of 
information is characterized by the change of the distin- 
guishability between a pair of quantum states which, in 
turn, is linked to the change of the trace distance between 
these states. We have also argued why the trace distance 
represents the most suitable distance measure for quan- 
tum states to be used in this context. Furthermore, since 
we have developed a genuine quantitative measure, the 
results presented here also allow to compare the degree of 
non-Markovianity of different types of physical systems. 

It has been demonstrated that a nonzero measure for 
non-Markovianity requires the dynamical map to be non- 
divisible, a property which is thus necessary for the 
presence of memory effects in the open system dynam- 
ics. It has also been shown that Markovian semigroups 
and time-dependent Markov processes are divisible and, 
hence, lead to A/"($) = 0. The examples discussed here 
illustrate how the measure can be calculated for a given 
open system dynamics and that a nonzero measure for 
non-Markovianity is linked to the emergence of negative 
decay rates in the corresponding master equation. 

Our measure for non-Markovianity has a clear opera- 
tional meaning based on the interpretation of the trace 
distance in terms of the distinguishability of states, and 
suggests various ways to experimentally decide whether 
a system under study is non-Markovian. The measure- 
ment scheme discussed here has the great advantage that 
it does not presuppose any knowledge about the structure 
of the environment or about the system-environment in- 
teraction and, therefore, also gives valuable information 
on the theoretical modelling of the open system dynam- 
ics. If, for example, a substantial increase of the trace 
distance is observed experimentally, a mathematical de- 
scription of the dynamics through any equation describ- 
ing a Markovian or time-dependent Markovian process is 
excluded. This shows that our measure is a useful tool for 
the characterization of non-Markovianity, both in exper- 
iments on open systems and in their theoretical analysis 
and modelling. 

We have argued that the characteristics of the infor- 
mation exchange between the system and its environ- 
ment determine the degree of non-Markovian behavior in 
an open system. This exchange of quantum information 
has been defined here in very general terms through the 
change of the distinguishability of quantum states, and 
does not presuppose anything about the specific physi- 
cal carriers of the information, e.g., energy or particles. 
Moreover, the measure does not depend on any specific 
representation of the open system's dynamics. It there- 
fore opens the possibility to compare and assess differ- 
ent mathematical formulations of dynamical processes 
in their ability to describe memory effects, in order to 
understand better the mathematical description of non- 



Markovian quantum dynamics. 
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Appendix 

Here we present some details of the A-model studied 
in Sec. |III B| The weak-coupling master equation for this 
model is given by 



where 



dt 



p(t) = -iX 1 (t)[\a)(a\,p(t)}~iX 2 (t)[\a)(a\,p(t)} 



+71 (*) 
+72(0 



|6) (a\p(t)\a) (b\--{p(t),\a) (a\] 
\c) (a\p(t)\a) ( C \-Up(t),\a) (a\} 



where 



A,-(t) 
7i(t) 



t pOO 

ds / du J(u>) sin [(lu — Wj)s], 
Jo 

t poo 

ds / dojJ(u>) cos [(oj — u>i)s\. 
o Jo 



Introducing the definitions 
/(*) 

9i(jt) = / t d S7 ,( S ) e - [Dl(s)+D2(s)] , 



e -[D 1 (t)+D 2 (t)]/2 e -i[L 1 (t)+L 2 {t)] 



Di(t) 



dsji(s), Li(t) 



dsXi(s), 



the solution of the master equation can be represented 
as follows, 



Paait) 


= \f{t)\ 2 Paa(0), 




Pbb{t) 


= gx(t)p aa (0) + 


Pbb(0), 


Pcc{t) 


= g2{t)Paa{0) + 


Pcc(0), 


Pab{t) 


= f(t)Pab(0), 




Pac{t) 


= /(t)Pac(O), 




Pbc(t) 


- Pfcc(0). 





These equations define the dynamical map $(t, 0) of the 
A-model. Employing the results of Choi [25] one can 
check that a necessary and sufficient condition for the 
complete positivity of this map is given by 



ffi(i)>0, g 2 (t)>0. 



These conditions are satisfied for the parameters used in 
the simulations of Sec. IHI B[ 
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